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Abstract 

Clonally derived bacterial populations exhibit significant genotypic and phenotypic diversity that contribute to fitness in 
rapidly changing environments. Here, we show that serial passage of Salmonella enterica serovar Typhimurium LT2 (SfLT2) in 
broth, or within a mouse host, results in selection of an evolved population that inhibits the growth of ancestral cells by 
direct contact. Cells within each evolved population gain the ability to express and deploy a cryptic "orphan" toxin encoded 
within the rearrangement hotspot (rhs) locus. The Rhs orphan toxin is encoded by a gene fragment located downstream of 
the "main" rhs gene in the ancestral strain SfLT2. The Rhs orphan coding sequence is linked to an immunity gene, which 
encodes an immunity protein that specifically blocks Rhs orphan toxin activity. Expression of the Rhs orphan immunity 
protein protects ancestral cells from the evolved lineages, indicating that orphan toxin activity is responsible for the 
observed growth inhibition. Because the Rhs orphan toxin is encoded by a fragmented reading frame, it lacks translation 
initiation and protein export signals. We provide evidence that evolved cells undergo recombination between the main rhs 
gene and the rhs orphan toxin gene fragment, yielding a fusion that enables expression and delivery of the orphan toxin. In 
this manner, rhs locus rearrangement provides a selective advantage to a subpopulation of cells. These observations 
suggest that rhs genes play important roles in intra-species competition and bacterial evolution. 
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Introduction 

Bacteria often reside in complex communities such as biofilms 
in which cells from multiple species touch one another in a three- 
dimensional network [1]. These environments provide opportuni- 
ties for cellular interactions, yet the mechanisms underlying 
contact-dependent competition and cooperation have been largely 
unexplored until recently. A diverse family of YD-peptide repeat 
proteins mediates at least two distinct forms of contact-dependent 
competition in Gram-negative and -positive bacteria [2] . The Rhs 
(rearrangement hotspot) proteins of Gram-negative enterobacteria 
[3,4] are large (~ 1,400-1,700 residues) toxic effectors that appear 
to be exported through the type VI secretion machinery. Related 
WapA (wall-associated protein A) proteins from Gram-positive 
bacteria are somewhat larger (~2, 200-3, 600 residues) [5] and are 
likely exported through the general secretory pathway [2]. Rhs 
and WapA proteins are both characterized by sequence-diverse C- 
terminal regions (Rhs-CT and WapA-CT) that vary considerably 
between different strains of the same species. Analysis of several 
Rhs-CTs and WapA-CTs from Dkkeja dadantii 3937 and Bacillus 
subtilis subspecies revealed that these domains contain the toxin 
activities responsible for intercellular growth inhibition. All rhs and 
wapA genes are closely linked to small downstream open reading 



frames that encode Rhsl and Wapl immunity proteins, respec- 
tively. These immunity proteins are also sequence-diverse and only 
protect against their cognate Rhs-CT (or WapA-CT) toxins. Thus, 
Rhs and WapA represent related, yet distinct, delivery platforms 
for polymorphic toxin domains [2]. Because different strains 
typically express unique rhs-CT / rhsl (wapA-CT/ wapl) alleles, these 
systems collectively form a complex network of toxin/immunity 
pairs that are thought to mediate inter-strain competition for 
environmental resources [2]. 

The rhs loci of Enterobacteriacae often contain one or more 
additional rhs-CT/ rhsl gene pairs located downstream of the main 
rhs /rhsl pair. These modules have been termed "orphan" toxin/ 
immunity pairs, because the rhs-CT coding sequences resemble 
displaced fragments from full-length rhs genes [6]. Orphan rhs-CT 
genes often contain some coding sequence for portions of the 
conserved N-terminal regions, but orphan fragments are much 
smaller than full rhs genes and usually lack translation initiation 
signals. Therefore, it is unclear whether orphan rhs-CT genes are 
expressed, raising the question of whether these auxiliary elements 
are functional. Here, we show that repeated passage of Salmonella 
enterica serovar Typhimurium LT2 (SLT2) produces "evolved" 
lineages that deploy the orphan Rhs-CT toxin to inhibit the 
growth of ancestral cells. We provide evidence that the rhs locus 
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Author Summary 

Salmonella Typhimurium is a bacterium that causes 
intestinal diseases in a number of animals including 
humans. In mice, this pathogen invades tissues, causing 
symptoms similar to typhoid fever. In an effort to 
understand the evolution of this pathogen, we grew S. 
Typhimurium in either liquid broth or in mice for many 
generations and examined the resulting "evolved" strains 
to determine if they were different from the original 
"parent" culture. We found that many of these evolved 
strains inhibited the growth of the parent after they were 
mixed together, and that this growth inhibition requires 
that the evolved and parental cells are in close contact. 
Genetic analysis showed that this contact-dependent 
growth inhibition requires Rhs protein, which has a toxic 
tip. Salmonella is normally resistant to its Rhs toxin 
because it also produces an immunity protein that blocks 
toxin activity. However, evolved cells have undergone a 
DNA rearrangement that allows them to express a 
different Rhs toxic tip that inhibits growth of the parental 
cells, which lack immunity to it. This allows the evolved 
cells to outgrow the original parental cells. Our work 
indicates that populations of Salmonella are dynamic, with 
individuals battling with each other for dominance. 



undergoes rearrangement to fuse the rhs"""" and rhs-CT'^ " genes, 
thereby providing a mechanism to express and export the Rhs- 
CT orphan toxin domain. These results indicate that rhs rearrange- 
ment provides a selective advantage to a subpopulation of cells, 
suggesting that rhs plays an important role in clonal selection and 
bacterial evolution. 

Results 

In an effort to isolate SLT2 strains with increased fitness, we 
serially passaged cells for ~ 1 ,000 generations in LB medium [7] . 
Analysis of six independendy evolved cultures revealed that each 
lineage outcompeted ancestral &LT2 cells in co-culture experi- 
ments (Figures 1A & SI A). Remarkably, we observed the same 
competitive advantage in four of eight 6YLT2 lineages that were 
obtained by passage through multiple mouse hosts [8] (Figures 1A 
& SIB). This competitive advantage was not due to faster growth 
rate, because four of the evolved lineages grew more slowly than 
the ancestral strain (Figure S2). To further explore this phenotype, 
we tested whether evolved lineages inhibit ancestral cells in a 
contact-dependent manner. We co-cultured evolved and ancestral 
cells using trans-well culture dishes, in which the two populations 
are separated by membranes of different porosities [9]. The 
growth of ancestral cells was inhibited when the populations were 
separated by a cell-permeable 8.0 Urn filter, but not when cell 
contact was prevented with a 0.4 um filter (Figure IB). These 
results indicate that evolved cells must be in close proximity to 
target cells in order to inhibit growth. This phenomenon is 
reminiscent of Rhs-mediated growth inhibition, which we recently 
characterized for D. dadantii 3937 [2]. 5<LT2 contains a single rhs 
locus, which contains a full-length "main" rhs gene (STM0291) 
and an "orphan" rhs gene fragment (STM0292) (Figure 2). Both 
rhs genes are closely linked to small open reading frames 
representing potential rhsl immunity genes (Figure 2), although 
the predicted rhsF""" immunity gene found downstream of rhs™'" 1 is 
not annotated in the genome sequence NC_003197. To determine 
if the rhs region is responsible for the observed growth inhibition, 
we tested whether over-expression of either rhsr am or rhsl" f '" 1 " 



immunity genes provided protection against evolved &LT2 
lineages. Parental SLT2 cells overexpressing rhsT"'" were still 
inhibited by the evolved lineages, but overexpression of the 
rhsf^ gene fully protected targets from growth inhibition 
(Figure 1A). These data strongly suggest that evolved SLT2 cells 
gained the ability to deliver Rhs-CT orpllan toxin into neighboring 
cells. 

We next tested each rhsl rhsl gene pair to confirm that they 
encode functional toxin and immunity proteins. Nucleotides 
3608 to 4095 of rhs"" 11 " and nucleotides 269 to 741 of rhs-CT' 1 ''"'" 
were cloned under the control of the arabinose-inducible Pbad 
promoter. The predicted rhsl immunity genes were cloned using 
a compatible plasmid under control of the IPTG-inducible P ire 
promoter. These plasmids were then introduced into StUT2 
cells to evaluate toxin and immunity functions. Induction of 
either rhs-CT""" or rhs-CT' 1 ' 1 "'" in S(LT2 resulted in rapid 
growth arrest (Figure 3A). In each instance, growth inhibition 
was neutralized by expression of the cognate rhsl immunity 
gene. However, co-expression of non-cognate immunity genes 
did not alleviate growth arrest (Figure 3A), demonstrating that 
Rhsl mam and Rhsl" rphan immunity proteins are specific for their 
cognate toxins. We obtained essentially identical results upon 
expressing the rhs main toxin and immunity genes in E. coli cells 
(Figure 3B). These results indicate that Rhs-CT" rphan is capable 
of inhibiting bacterial growth and support a model in which 
evolved St\JT2 lineages deploy the orphan toxin to inhibit the 
ancestral strain. 

The rhs-CT* 1 ' 1 "'" sequence does not encode a full-length Rhs 
protein, raising the question of how this toxin is synthesized and 
exported from evolved cells. The rhs"""" and rhs-CT" 1 ' 1 ""' coding 
regions share 95% sequence identity over 522 base-pairs (Figure 2), 
raising the possibility that homologous recombination in the 
evolved lines generates a new full-length rhs gene that encodes the 
Rhs-CT orpllan toxin domain [10]. Bacteria expressing this Rhs 
chimera would have a growth advantage if rhsT rt,lmn expression is 
low in ancestral cells. However, the proposed recombination event 
would also delete the rhsT m " gene, rendering the evolved cells 
sensitive to inhibition by siblings expressing the main Rhs-CT 
toxin. Therefore, we hypothesized that rhs recombination occurs 
subsequent to duplication of the locus such that evolved cells retain 
the rhsT""" immunity gene (Figure 2). To test this hypothesis, we 
analyzed chromosomal DNA from evolved and ancestral lineages 
by Southern blot. DNA was digested with Hindi, which cleaves 
between the rhsl 1 """ and rhs-CT"^" coding sequences, and probed 
with a labeled DNA fragment that specifically hybridizes to rhs"""" 
(Figure 2). We detected a unique junction fragment representing 
fusion of rhs-CT^ 1 " 1 " to the upstream rhs"""" gene in StLT2 lineage 
2, which displayed the highest level of growth inhibition of all 
lineages (Figures 1 & 4A). The wild-type rhs locus was also detected 
in lineage 2 (Figure 4A), which is consistent with rhs region 
amplification, but may also indicate distinct populations of 
recombinant and non-recombinant cells. Orphan rhs recombi- 
nants were not detected in the other evolved lineages by Southern 
blot analysis (Figure 4A). Because the growth inhibition phenotype 
varied in magnitude between the different evolved strains, it is 
possible that only a fraction of the evolved StLT2 cells are rhs 
recombinants. If so, then the proportion of recombined rhs loci in 
the DNA sample may be below the detection limit of Southern 
analysis. Therefore we analyzed each evolved lineage with 
quantitative real-time PCR (qPCR) to measure the relative levels 
of rhs ma, "-rhs°' l ' hm ' junction sequences. All five of the evolved lineages 
contained 10- to 1,000-fold more rhs"""" -rhs or,,ha " junction than 
ancestral $LT2 (Figure 4B), consistent with the ability of these 
strains to deploy Rhs-CT orphan toxin. 
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Figure 1. Evolved SA.T2 cells inhibit the growth of ancestral cells. A) The indicated evolved Sfl_T2 lineages were co-cultured with the 
ancestral strain for 24 h in broth. Viable cell counts for each population were determined as colony forming units and these data were used to 
calculate the competitive index as described in Methods. Each evolved lineage was competed against ancestral cells (black bars), ancestral cells 
overexpressing rhsl mam (light grey bars) and ancestral cells overexpressing rhsf" phan (white bars). Reported values represent the mean ± SEM for at 
least three independent experiments. B) The growth inhibition activity of evolved 5fl_T2 requires cell-cell contact. Evolved cells were co-cultured with 
ancestral cells in adjacent wells of a frans-well incubation chamber. Culture chambers were separated by membranes containing 0.4 |im or 8 |j.m 
pores as indicated. Reported competitive indices represent the mean ± SEM for three independent experiments. 
doi:10.1371/journal.pgen.1004255.g001 
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Figure 2. Proposed mechanism of rhs locus rearrangement in evolved inhibitor cells. The 5fLT2 rhs locus contains a full-length rhs mam 
gene (STM0291) and an rhs-CT° rphan gene fragment (STM0292) together with associated rhsl immunity genes. Duplication of the rhs locus would 
provide the opportunity for subsequent homologous recombination between the 522 bp region of near sequence identity (95%) between rhs maln 
and rhs-CI° rphan (depicted as diagonal hatched regions). The binding site of the Southern blot hybridization probe is indicated by the double ochre 
bars. Primer binding sites for PCR amplification of the rhs""" n /rhs-CTorphan junction are indicated by convergent horizontal arrows. 
doi:1 0.1 371 /journal.pgen.1 004255.g002 
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Figure 3. The 5A.T2 rhs locus encodes cognate toxin/immunity pairs. A) Expression of rhs-C~T"" n or rhs-CJ° rphan inhibits the growth of SfLT2 
cells. Expression the plasmid-borne rhs-CT genes was induced by addition of L-arabinose at the indicated time, and cell growth was monitored by 



measuring the optical density at 600 nm (OD 60 o). The cells also co-expressed either rhsl mal " (dark squares) or rhsf 



prphan 



mmunity genes (light grey 



diamonds) from IPTG-inducible promoters. Growth is compared to control cells that carry the empty vector plasmids (triangles). B) Expression of rhs- 
Cr°'" n or rhs-cr rpha " inhibits the growth of E.coli cells. The rhs-CT and rhsl ge nes were expressed in E. coli cells from the same plasmids described in 
panel A, and growth monitored by measuring the OD 60 o of the cultures. 
doi:1 0.1 371 /journal.pgen.1 004255.g003 



Because only a fraction of the passaged cells appeared to display 
growth inhibitory activity, we asked whether inhibitor-cell clones 
could be isolated from each population. As a control, we first 
isolated colonies from an overnight culture of the ancestral strain 
and tested these clones for growth inhibition activity. None of the 
ten ancestral clones tested were inhibitory, suggesting that the 
proposed rhs rearrangements occur at low frequency. By contrast, 
approximately 30-90% of the clones isolated from the culture- 
evolved lineages and ~20% of the clones from mouse-evolved 
lineage 1 showed inhibition activity against ancestral cells 
(Figure 5A). However, no inhibitor clones were isolated from 
mouse-evolved lineage 2 (Figure 5A). Strikingly, the inhibition 
activity of these clones varied considerably. For example, 
competitive index values ranged from 10 to 10 J for 
competitions between ancestral cells and inhibitory clones isolated 
from evolved lineage 2 (Figure S3). Although their potencies 
varied, it appears that each inhibitor-cell clone deployed the Rhs- 
Qrporphan tox - n because ancestral cells could be protected through 
over-expression of rhsf' 1 ' 1 "'", but not rhsl" 1 "" (Figure 5B). The 
presence of DNA fragments corresponding to both ancestral and 
recombinant rhs loci in lineage 2 (Figure 4A) suggests that either 



the rhs region was duplicated or there are distinct populations of 
recombinant and non-recombinant cells. In the latter case, single 
colonies isolated from the inhibitory lineages would contain only 
the rhs- r/z.s-CT orph ' m junction and not the rfe-CT™™ sequence. 
However, PCR analysis of the single colonies with inhibitory 
activity in Figure 5B showed that each contained both ancestral 
and recombinant rhs loci. In addition, sequence analysis of the 
recombinant PCR product verified that recombination occurred 
between the regions of homology shared by rhs"""" and rhs-CT"^"" . 
Together, these data demonstrate that the evolved populations are 
heterogeneous with respect to Rhs-CT orphim mediated inhibition 
activity. Furthermore, these results suggest that the inhibition 
phenotype of a given culture may be due entirely to a minor 
subpopulation of potent inhibitor cells. 

Because inhibitor cells represent a subpopulation in the evolved 
cultures, the other non-recombinant cells in the cohort are 
presumably resistant to the Rhs-CT orph ' m toxin. To test this 
hypothesis, we isolated non-inhibitory clones from each of the 
evolved cultures and tested them in competition co-cultures 
against their respective evolved lineages. As predicted, each of the 
non-inhibitory clones was either fully- or partially-resistant to its 
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Figure 4. Evidence for rhs rearrangement in evolved 5/LT2. A) Southern blot analysis of evolved lineages. Arrows indicate positions of the 
2.8 kbp Hindi restriction fragment containing the ancestral rhs region and the 4.7 kbp restriction fragment resulting from recombination between 
rhs maln and rhs-Cr" phan (see Figure 2 for model). B) Real-time qPCR analysis of rhs lr "" n /rhs-Cr rphan - cr recombination junctions in evolved 5fLT2. The 
levels of rhs""""/rhs-CT° rpha " junction products are expressed relative to amplified products of a control locus (bamA). Positive control cells are 
engineered Sfl_T2 that contain a chromosomal deletion fusing rhs"""" to rhs-Cl° rphan '. 
doi:1 0.1 371 /journal. pgen.1004255.g004 



cohort lineage (Figure 5C). These cells likely carry uncharacterized 
resistance mutations that may prevent cell-cell contact, block the 
delivery of Rhs-CT orphan toxin,, or increase the immunity of these 
cells to Rhs-CT toxin. 

To directly detect Rhs-CT orphan expression in the evolved 
lineages, we examined cells by immunofluorescence microscopy 
using polyclonal antibodies against the Rhs-CT orphan toxin. Rhs- 
Qrporphan ant ;g en was detected on the surface of some cells within 
evolved lineages 1, 2 and 3 as well as mouse-evolved lineages 1 
and 2 (Figures 6A & S4). In contrast, the Rhs-CT" rphan signal was 
undetectable on the surface of both ancestral SLT2 cells and cells 
carrying a deletion of the rhs-CT"^ " (Figures 6A & S4). We then 
quantified the fraction of cells with Rhs-CT orphan antigen on the 
cell-surface using flow-cytometry. Evolved lineages showed a 2- to 
20-fold increase in the fraction of Rhs-CT orpllan -positive cells 
compared to ancestral SLT2 cells (Figures 6B & S5). Mouse- 
evolved SLT2 showed very low expression of Rhs-CT orphan 
antigen on cell surfaces (Figure 6B), consistent with the modest 
growth inhibition observed for these lineages (Figure 1). Based on 
Southern blot and RT-qPCR analyses, it seems likely that surface 
expression of Rhs-CT orphan requires rhs locus rearrangement to 
generate a chimeric rhs gene. In accord with this conclusion, we 
also found that over-expression of rhs-C'I° ,t ' > " m from a multicopy 
plasmid does not increase Rhs-CT orpllan antigen levels on the cell 
surface (Figures 6B & S5). Therefore, we sought to detect the 
predicted Rhs fusion protein using antisera to the Rhs-CT orphan 
toxin. Western blot analysis revealed an immuno-reactive protein 
at — 150 kDa in culture evolved lineages 2 and 3 (Figure 6C). 
This product corresponds to the expected size of the Rhs fusion 
protein. Moreover, we were unable to detect the 29 kDa product 
encoded by rhs-CT"® in the ancestral and evolved lineages 
(Figure 6C). Together, these data strongly suggest that the rhs- 
Qjerphan reac y n g f mme must recombine with rhs"" 11 " to be 
expressed. 

Discussion 

The results presented here show that serial passage of SLT2, in 
either laboratory media or within a natural host, leads to 
enrichment of cells that express Rhs-CT orphan toxin. Analysis of 



the rhs locus indicates that evolved cells undergo recombination 
between rhs"""" and rhs-CT"^"' ', forming a gene fusion that allows 
the Rhs-CT orphan toxin domain to be deployed. Rearranged rhs 
genes are detected at low levels within the evolved populations, 
indicating that only a fraction of cells are recombinant inhibitors. 
A number of observations argue that this subpopulation of cells is 
responsible for growth inhibition activity. First, the relative 
competitive advantage of each evolved lineage is correlated with 
its level of recombinant rhs junctions and surface expression of 
Rhs-CT orpllan antigen. More importantly, ancestral cells are fully 
protected when they over-express the rhsl" t ' h '"' immunity gene. 
Because Rhs immunity proteins are highly specific for their 
cognate toxins, this latter result demonstrates that Rhs-CT orphan 
toxin is indeed deployed by the evolved lineages. This result also 
indicates that ancestral SLT2 cells do not normally express 
rhsF® immunity genes under laboratory conditions. The number 
of inhibitor cells within each lineage is not known, but can be 
estimated to be <2% of the population based on flow cytometry 
measurements of Rhs-CT orphan antigen on cell surfaces. However, 
we note that this assay may underestimate the actual number of 
recombinant inhibitor cells because Rhs effectors are likely 
exported through type VI secretion systems [2,11]. Although 
recent studies indicate that the N-terminal PAAR domain found 
within many Rhs proteins forms the tip of the type VI injection 
structure [12], other structural studies show that Rhs-peptide 
repeats form a chamber capable of encapsulating toxin domains 
[13]. Therefore, much of the Rhs-CT olphan antigen may be 
inaccessible to antibody until it is delivered to target cells. In 
accord with this model, we only detect Rhs-CT orphan where two 
bacteria make contact with one another and never on the surface 
of individual cells. Regardless of the absolute number of 
recombinants or rhs expression levels, our results suggest that a 
small number of inhibitor cells are capable of inhibiting a large 
excess of ancestral cells. The same phenomenon has been 
observed during bacterial contact-dependent growth inhibition 
(CDI), in which each CDI + ceU is able to inhibit 100-1,000 target 
cells over a few hours [9]. Presumably, the unstructured 
environment in shaking broth culture promotes a series of 
transient cell-cell interactions, thereby enabling toxin delivery to 
multiple ancestral cells. 
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Figure 5. A subpopulation of evolved cells has growth inhibition activity. A) Two sets of independent clones were isolated twice from the 
evolved lineages and tested for growth inhibition activity against ancestral cells. The percentage of evolved clones with inhibition activity is shown. 
Reported values represent the mean ± SEM for at least two independent experiments. B) Growth inhibition activity of isolated evolved clones. Clones 
from each evolved lineage were competed against ancestral cells (black bars), ancestral cells overexpressing rhsf" am (light grey bars) and ancestral 
cells overexpressing rnsP' pha " (white bars). Reported values represent the mean ± SEM for at least three independent experiments. C) Growth 
inhibition activity of evolved lineages towards non-inhibitory clones. Evolved lineages were co-cultured with ancestral cells (black bars) and two non- 
inhibitory clones isolated from the evolved cultures (light grey and white bars). Competitive indices represent the mean ± SEM for at least three 
independent experiments. 
doi:1 0.1 371 /journal.pgen.1 004255.g005 



Chromosomal duplications and amplifications occur frequently 
in bacteria, typically at rates of about 0. 1 % per generation for any 
given locus [14,15]. However, there is a cost to maintaining 
amplified regions, and gene duplications are lost during segrega- 
tion at frequencies up to 10% per generation [16,17]. Therefore, 
positive selection is required to retain multiple gene copies. If the 
amplified region can be stabilized, then the additional gene copy 
can diverge towards a new function, thus providing a mechanism 
for evolution [16,18]. Rearrangement of rhs loci represents a 
previously unrecognized mechanism for bacteria to exploit 
chromosomal amplifications for adaptation. We propose that, 
subsequent to duplication, homologous recombination occurs 
between rhs"""" and rhs-C7 1>r> ' > "'" to generate a novel chimeric rhs 
element. This recombination would necessarily delete one copy of 
the rhsF""", but the other copy would remain and ensure that 
recombinant cells retain immunity to the Rhs m "* m toxin should it 
be deployed by neighboring non-recombinant siblings. This model 
also predicts that evolved recombinant cells could undergo 



homologous recombination to restore the original rhs locus (see 
Figure 2, reverse of the duplication step). Thus, rhs rearrangement 
could be exploited transiently under conditions where it confers a 
selective advantage, but rapidly revert back to the ancestral 
genotype as environmental circumstances dictate. 

Analysis of over 150 Salmonella genomes shows that rhs-CT tasasx 
sequences are diverse with at least 57 distinct sequence types 
(Figure S6A & Table SI). This is a common feature of rhs genes in 
other bacteria as well and suggests that Rhs mediates inter-strain 
competition. All Salmonella serovars contain at least one rhs gene, 
located on pathogenicity islands SPI-6 or SPI-19 [19,20]. 
Approximately 50% of these serovars contain at least one 
predicted rhs orphan sequence, with some strains containing as 
many as eleven modules. There is generally high conservation of 
Rhs-CT main and Rhs-CT orphim . sequences within a given serotype. 
For example, all sequenced Typhi isolates contain the same Rhs- 
CT mam and Rhs-CT orphan . sequences, whereas these CT sequence 
types are only found in one other serotype, thus suggesting that 
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Figure 6. Expression of Rhs-CT orphan in evolved cells. A) Immunofluorescence analysis of ancestral 5fl_T2 and culture-evolved lineage 3 with 
antibodies against Rhs-CT orphan . Scale is 10 u.mx10 am for each image. B) Quantification of surface-expressed Rhs-CT° rphan . Cells were labeled with 
Rhs-CT° rphan antisera and analyzed by flow cytometry as described in Methods. Reported values represent the mean ± SEM for three independent 
experiments with 50,000 events recorded per sample. C) Immunoblot analysis of Rhs-CT orphan . Proteins were isolated from ancestral and evolved cells 
and analyzed by immunoblot using antibodies against Rhs-CT° rphan . Regions corresponding to predicted migration positions of Rhs-CT orphan 
(—31 kDa) and the chimeric fusion protein (—151 kDa) are shown. 
doi:1 0.1 371 /journal.pgen.1 004255.g006 



different toxins are linked to serotype and/ or the type of infection. 
However, orphan rfc-CTsequences in one serotype can be present 
within the main rhs gene of another serotype. For example, the 
&LT2 Rhs-CT orphan toxin studied in this work is part of the full- 
length main Rhs in Salmonella enterica serovar Saintpaul SARA23 
and some Newport isolates (Figures S6A & S6B). These 
observations and the association with horizontally transferred 
elements suggest that rhs genes are exchanged between different 
serovars and contribute to the evolution of toxin diversity. 

Given that Rhs toxins are encoded on pathogenicity islands, it 
seems likely that these systems also play important roles in 
Salmonella growth and fitness during pathogenesis. Indeed, S/LT2 
mutants lacking a chromosomal region containing rhs-CT"^ 1 " 1 " are 
outcompeted by wild-type cells in mice [21], and SSL 1344 
mutants lacking rhs are completely attenuated in pig and cattie 
models of infection [22]. These observations raise the possibility 
that rhs locus rearrangement occurs commonly during infections. 
Intriguingly, SthT2 produces distinct intracellular infection foci, 
each originating from one or only a few clones [23]. Similarly, 
analysis of mice orally infected with Yersinia pseudotuberculosis 
indicates that only a few bacterial clones are able to disseminate 
from the intestines to the spleen and liver [24] . Clonal invasion has 
also been reported for Yersinia enterocolitica infections [25], but the 
mechanisms underlying these apparent dissemination bottlenecks 
are unknown. Most Yersinia species contain rhs loci with associated 
orphan gene pairs, raising the possibility that clonal expansion 
through rhs recombination and growth selection may be a general 
feature of many enterobacterial infections. Rearrangement could 
function as a stochastic switch that enables some cells to deploy 
Rhs-CT" rph ' ln and thereby "differentiate" into cells that are 
specialized for tissue invasion or immune modulation. Although 
Rhs-mediated inhibition clearly occurs between bacteria, it is 
also possible that Rhs toxins act directly as virulence factors. The 



C-terminal region of RhsT from Pseudomonas aeruginosa was recendy 
shown to be delivered into mouse host cells [26]. In the process, 
the Rhs fragment activates the inflammasome and contributes to 
pathogenicity. 

Materials and Methods 

Strains and growth conditions 

Bacterial strains were derived from Salmonella enterka serovar 
Typhimurium LT2 (iS(LT2) and are listed in Table S2. Bacteria 
were grown in LB medium [27] supplemented with 50 mM 
potassium phosphate (pH 7.3). Bacteria were incubated at 37°C 
with shaking at 200 rpm. Where appropriate, media were 
supplemented with antibiotics at the following concentrations: 
ampicillin (Amp), 200 mg/L; chloramphenicol (Cam), 17 mg/L; 
kanamycin (Kan), 80 mg/L; and tetracycline (Tet), 5 mg/L. Six 
independent lineages of SLT2 were obtained by serial passage for 
~ 1 ,000 generations [7] . Each lineage was passaged daily by 
dilution of 1.5 uL of overnight culture into 1.5 mL of fresh LB 
medium. Each evolved lineage was sampled periodically (100-150 
generations) and stored at — 80°C. 

All growth competitions were conducted using ancestral StUI2 
marked with the jlhCy.cat allele, which confers Cam resistance. 
Non-inhibitory clones isolated from the evolved cultures were 
transduced with the JlhCr.cat allele prior to testing for resistance. 
Ancestral and evolved cells were co-cultured in LB medium 
supplemented with 50 mM potassium phosphate (pH 7.3) at 37°C 
with shaking. At time Oh, ~ 1 0 6 cfu ( 1 |jL of overnight culture) 
from evolved and ancestral cultures were suspended in 2 mL of 
fresh LB (pH ~7.3) and plated for viable cell counts before shaking 
incubation for 24 h at 37°C. After 24 h of co-culture, viable cell 
counts were determined by plating onto LB agar (to enumerate 
evolved and ancestral cells) and LB agar supplemented with Cam 
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(to enumerate ancestral cells). The competitive index was 
calculated as the ratio of ancestral: evolved cells at time 24 h 
divided by the cell ratio at 0 h. Ancestral SLT2 JlhCr.cat cells were 
also supplemented with either rhsl" am or rhsl" 1 ' 11 "" on the 
chromosome under control of the lac promoter and on plasmid 
pBR322 under the kt promoter. Chromosomal rfo/ >rphan and 
plasmid-borne rhsl" plmn individually provided partial protection 
against the evolved lineages (data not shown), but both copies were 
required for full immunity. Proximity-dependence of growth 
inhibition was determined as described previously [9]. Cells were 
grown to ODgoo ~0.3, then transferred to a trans-well culture plate 
(BD diagnostics) that separates the two populations with filter 
containing 0.4 urn (no-contact) or 8.0 u,m (contact) pores. Trans- 
well culture plates were seeded at an evolved:ancestral cell ratio of 
1:1 and incubated at 37°C with shaking for 24 h. Cultures were 
then plated onto selective media to determine viable cell counts 
and to calculate competitive indices. 

Construction of plasmids and chromosomal inserts 

All oligonucleotides used in this study are presented in Table S3. 
The rhsF""" and rhsl" 1 ' 1 '"" genes were amplified from ancestral 
SLT2 chromosomal DNA using oligonucleotides 2337/2338 and 
2340/2544 (respectively) and ligated to plasmid pBR322 using 
EcoRV and Sail restriction sites. The immunity genes were also 
placed under the lac promoter at the glmS locus using bacterio- 
phage X Red-mediated recombination [28] . Integration constructs 
containing rhsl genes flanked by a Kan-resistance cassette and 
glmS-derived homology regions were constructed by overlapping 
end-PCR as described previously [29]. The following primer pairs 
were used to amplify: upstream glmS homology (2666/2676), lac 
promoter (2677/2678 for rhsT m and 2677/2682 for rhsl" pha \ 
rhsF' 1 " 1 (2679/2680) or rhsl" 1 ' 1 "" 1 (2683/2684), Kan-resistance 
cassette (2618/2619) and downstream glmS homology (2681/ 
2667). The final PGR product was electroporated into SLT2 cells 
that express Red recombinase proteins, and transformants were 
selected on LB supplemented with Kan. Integrated immunity 
genes were verified by PGR analysis using primers 2666/2667 and 
subsequent DNA sequencing. The JlhCr.cat and STM0292::fcm 
alleles were generated by PCR using primers 2436/2437 and 
2410/2490 to amplify the cat/kan cassettes of plasmids pKD3 and 
pKD4, respectively. Each PCR product was integrated into the 
SLT2 chromosome by Red-mediated recombination. 

To evaluate toxin activity and the specificity of immunity, 
individual rhs-CT and rhsl sequences were cloned under the control 
of inducible promoters on compatible plasmids. The rhs-CT mam 
and rhs-CT' 1 ' 1 ""' coding sequences were amplified with primers Sty- 
rhs(E1203)-Nco/Sty-rhs-Xho and Sty-rhs(E1203)-Nco/Sty-orph- 
rhs-Xho (respectively) and ligated to plasmid pCH450 [30] using 
Ncol and Xhol restriction sites. The rhsl™" 1 and rhsF* hm genes 
were amplified and ligated to a derivative of plasmid pTrc99A 
using Kpnl and Xhol restriction sites. Rhs-CT orpllan was expressed 
and purified as a non-toxic variant fused to His 6 -tagged 
thioredoxin. The hisg-trxA sequence was amplified from plasmid 
pSH21P::fr^4 [31] using primers pET-Sph and trxA-Bam-TEV- 
Kpn. The product was digested with Sphl/BamHI and ligated to 
plasmid pET21b to generate plasmid pSH21P::frx/TZEK The 
coding sequences for Rhs-CT orphan residues 112-247 and RhsI or " 
phan were amplified using primers Sty-rhs(D1225)-Kpn/Sty-orph- 
rhsI-Xho) and the His208Ala mutation made by mega-primer PCR 
using oligonucleotide Sty-CTol-H208A. The final product was 
digested with KpnI/XhoI and ligated to plasmid pSH21P::<rai4- 
TEV to generate plasmid pCH 10068. The resulting construct was 
used to overproduce His 6 -TrxA-Rhs-CT(H208A) orphin fusion 
protein. 



Chromosomal DNA analysis 

Chromosomal DNAs were isolated using the Sigma genomic 
DNA kit and digested with Hindi restriction endonuclease. 
Digested DNAs were resolved by electrophoresis on 0.7% agarose 
gels at 34V for 15 h and blotted onto nylon membranes by 
capillary transfer. A probe to nucleotides 2969-3 1 28 of rhs"""" was 
generated by PCR using oligos 2226/2227 and labeled with re- 
labeled using the Prime-It Random Primer Labeling Kit (Agilent 
Technologies). Southern blots were visualized by phosphor 
imaging. Fragment sizes were calculated using a standard curve 
based on Hindlll digested X ladder (New England Biolabs, USA) 
run on the same gel. The proportion of rhs recombination 
junctions was determined by quantitative real-time PCR (qPCR) 
using oligonucleotides 2226/2231 using the cycle threshold C t - 
value method according to the manufacturer (Bio-Rad). Fluores- 
cence was monitored on-line using the MyIQ_ iCycler real-time 
PCR system (Bio-Rad). The rhs-rhs mt ' hm junction DNA levels were 
calculated relative to bamA DNA (oligos 1981/1990) in each 
sample and normalized to the level of junction DNA in ancestral 
cells. 

Antiserum preparation and immunoblot analysis 

His 6 -TrxA-Rhs-CT(H208A) orphan fusion protein was overpro- 
duced in E. coli CH2016 and purified by Ni 2+ -affinity chroma- 
tography as described [32]. The Rhs-CT(H208A) orpllan domain 
was released by TEV protease digestion and used for antiserum 
production in rabbits (CoCalico Biologicals). Non-specific anti- 
bodies were removed by incubation with carbonyl-diimidazole- 
activated agarose beads linked to soluble protein from E. coli 
strain CH2016 [33]. Briefly, protein-linked beads were resus- 
pended in 0.5 mL of antiserum (1:5 dilution) and mixed by 
rotation for 1 h at room temperature followed by additional 
incubation for 3 h at 4°C. This process was repeated at least four 
times with fresh beads. 

Evolved lineages were grown to mid-log phase in LB medium 
supplemented with 50 mil potassium phosphate (pH 7.3) and 
cells were collected by centrifugation and frozen at — 80°C. Cell 
pellets were resuspended in NuPage LDS-sample buffer (Invitro- 
gen) at 70°C and treated with benzonase to degrade nucleic acids. 
Cell lysates were run on 3-7% NuPage Tris-acetate gradient gels 
(Novex) for the detection of the Rhs mam -Rhs-CT orphan chimera, or 
on 4-10% Precise Tris-glycine gradient gels (Thermo Scientific) to 
detect Rhs-CT orpllan . Gels were electrotransferred to nitrocellulose 
membranes and the blots incubated with polyclonal antisera 
against Rhs-CT olphan (1:1,000 dilution) and secondary anti-rabbit 
800CW antiserum (1:10,000 dilution). Immunoblots were visual- 
ized using an Odyssey CLx Infrared Imaging System (LI-COR). 

Immunofluorescence microscopy and flow cytometry 

Cells were incubated overnight with 4% formaldehyde in 
0.15 M phosphate buffered saline (PBS, pH = 7.2). Cells were 
washed three times with PBS and incubated with polyclonal 
antibodies to Rhs-CT orphan (1:50 dilution) for 30 min. Cells were 
washed with PBS before incubation with secondary anti-rabbit 
Alexa-Fluor 480 antibodies (1:500 dilution) (Invitrogen) for 30 min 
on ice. After washing with PBS, cells were applied to poly-D-lysine 
coated slides, treated with Fluoro-gel II/DAPI (Electron Micros- 
copy Sciences) and visualized by fluorescence microscopy. The 
fraction of evolved cells expressing Rhs-CT orphan on the surface 
was determined by flow cytometry. Antibody-labeled cells were 
analyzed (50,000 events for each sample) with an Accuri C6 flow 
cytometer with gates set to include bacteria-sized particles. StUT2 
L\rhs-CT ,phan cells were used to assess non-specific binding of 
the Rhs-CT <>rphan antisera. The fraction of cells with surface 
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Rhs-CT orphan antigen was calculated as the ratio of green 
fluorescent particles in the population after subtracting back- 
ground fluorescence observed with S(LT2 Arhs-CT" phan cells. 

Supporting Information 

Figure SI Evolved SiLT2 outcompete ancestral cells. The 
indicated evolved 5/LT2 lineages were co-cultured with the 
ancestral strain for 24 h in broth. Viable cell counts for each 
population were determined as colony forming units and these 
data were used to calculate the competitive index as described in 
Methods. A) Culture-evolved lineages after 1000-generations of 
growth in LB were competed against ancestral wild type cells. B) 
Mouse evolved lineages after 150-generations of growth in mice 
were competed against ancestral cells. Reported values represent 
the mean ± SEM for at least three independent experiments. 
(PDF) 

Figure S2 Growth rates of ancestral and evolved SILT2 strains. 
The growth rates of evolved lineages are expressed relative to the 
growth rate of ancestral cells, which was set to 1 . Reported values 
represent the mean ± SEM for at least three independent 
experiments. 
(PDF) 

Figure S3 Variability of growth inhibitory activities of evolved 
inhibitor clones. Culture-evolved lineages were streaked on LB 
agar plates to obtain individual colonies. Ten colonies from each 
lineage were competed against the ancestral SLT2 strain as 
described for Figure 1. Reported values represent the mean ± 
SEM for at least two independent experiments. The hatched lines 
in each panel indicate an arbitrary cut-off (C.I. = 10 ) for 
whether a clone was considered to express growth inhibitory 
activity or not. Clones with C.I. error bars that cross the hatched 
line were considered to express growth inhibitory activity. 
(PDF) 

Figure S4 Evolved cells express Rhs orphan ~ CT on the cell surface. 
Immunofluorescence analysis of ancestral cells and evolved 
lineages using antibodies against Rhs-CT orph ' ln . Non-permeabi- 
lized cells were fluorescendy labeled with antibodies to Rhs- 
CT „ q ,ha„ as described in Methods. The Arhs-CT rphax cells carry a 
deletion of the rhs-CT"'' m gene. P ara ::Rhs-CT olplian carry a 
plasmid encoded Rhs-CT orphan under an arabinose inducible 
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